The network topology of a potential energy landscape: A static scale-free network 
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Here we analyze the topology of the network formed by the minima and transition states on the 
potential energy landscape of small clusters. We find that this network has both a small-world and 
scale-free character. In contrast to other scale-free networks, where the topology results from the 
dynamics of the network growth, the potential energy landscape is a static entity. Therefore, a fun- 
damentally different organizing principle underlies this behaviour: The potential energy landscape 
is highly heterogeneous with the low-energy minima having large basins of attraction and acting as 
the highly-connected hubs in the network. 

PACS numbers: 89.75.Hc,61.46.-|-w,31.50.-x 



Energy landscapes have been at the forefront of many 
of the recent theoretical develnpnients in our understands 
ing of biomolcculesp clustersoa and the glass transition.cl 
For example, this research has. provided important new 
insights into how proteins foldU and the origin of the un- 
usual properties of supercooled liquids, such aa.lhe dis- 
tinction between "strong" and "fragile" liquids .BEl There 
has also been a surge q£ interest in modelling com- 
plex systems as networksja inspired by Watts and Stro- 
gatz's discovery that many networks behave as "small 
worlds" .u Intriguingly, a divetse range of|-such net- 
works, e.g. the-jAforld-wide W|eb,Ej the internet,lliJ scientific 
coUaboratianPl and citationll3 networks, and biochemical 
networks, Eitj also have a "scale-free" topology, where 
the distribution of the number of connections to each 
node, the degree, follows a power law. This topology re- 
sults from the dynamics of the network growth in these 
systemsJlj Here we draw these two strands of research to- 
gether by applying the techniques of network analysis to 
probe the global structure of potential energy landscapes 
of clusters. 

The potential energy landscape is a multi-dimensional 
surface representing the dependence of the potential en- 
ergy on the positions of all the atoms of the system. For 
a system with many atoms the landscape will have a 
complex topography with higher-dimensional analogues 
of mountain ranges, valleys and passes. Although the po- 
tential energy landscape determines the system's struc- 
ture, thermodynamics and dynamics, the nature of this 
relationship is complex. A particularly successful means 
of elucidating this relationship is tlie inherent structure 
approach of Stillinger and Weber|i3 in which the land- 
scape is divided into basins of attraction surrounding 
each minimum (See Fig. This partition provides a 
natural way to describe the dynamics of the system, be- 
cause, except at high temperature, the system spends 
most of the time vibrating in the well surrounding a min- 
imum and only occasionally hops to a different well by 
passing over a transition state. The interbasin dynamics 
can then be represented as a walk on a network whose 
nodes correspond to the minima and where edges link 
those minimum which are directly connected by a tran- 
sition state. Figure |l] provides an illustration of such an 




FIG. 1: (a) A model two-dimensional energy surface, (b) A 
contour plot of this surface illustrating the inherent structure 
partition of the configuration space into basins of attraction 
surrounding minima. The basin boundaries are represented 
by red lines, and the minima and transition states by blue and 
green points, respectively, (c) The resulting representation of 
the landscape as a network. 



inherent structure network (ISN) for a two-dimensional 
energy surface. 

The ISN should provide the starting point for an en- 
ergy landscape view of the global dynamics of a system. 
Indeed, as is increasingly being done, it is relatively easy 
to calculate the dynairu.cs . from this network using a mas- 
ter equation approach .cIBEj However, fundamental ques- 
tions about the topology of the ISN have received little 
attention. By contrast the global topography ofj-energy 
landscapes has been the focus of much researchBa'B To 
give one example, this emjjhasis has revealed that when a 
landscape is like a 'funnel'El the system is guided towarcLa 
the global minimum, be it, the native state of a protein,Ll 
an ordered nanoparticlc£3 or a bulk crystal. However, 
the topological aspects of this idea remain open despite 
their importance: if the average number of steps to reach 
the global minimum from an arbitrary starting minimum 
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FIG. 2: (a) The dependence of the average separation be- 
tween nodes (in steps) on the size of the network, A^min- (b) 
The size dependence of the clustering coefficient, c, compared 
to that for a random graph, where c is the fraction of the pairia 
of nodes with a common neighbour that are also connected.u 
The data points are for Lennard- Jones clusters with from 7 
to 14 atoms. 



scales unfavourably with size, the location of this struc- 
ture would become significantly hindered at large size. 

To characterize the topology of the ISN we study small 
clusters for which we are able to locate nearly all the 
minima aBd-|transition states on the potential energy 
landscapeflO The atoms of the cluster interact with a 
Lennard-Jones potential, which provides a reasonable de- 
scription for rare gas clusters. The numbers of minima 
and transition states are expected to increas e , rou ghly as 
Nmin ~ and Nts « iVe"^, respectively,BKa where 
N is the number of atoms in the cluster. Therefore, the 
largest network that we are able to consider is for a 14- 
atom cluster for which we have located 4196 minima and 
87 219 transition states 

Small-world networks have characteristics typical of 
both random graphs and lattices. The average separation 
between nodes scales logarithmically with network size, 
while the network is highly clustered, i.e. any two neigh- 
bours of a node are also likely to be connected. From 
Fig. 1^ it is clear that the ISNs for the clusters show both 
these features and so are small worlds. The clustering is 
unsurprising given that the connections between basins 
on a potential energy landscape are based on adjacency 
in configuration space,E3 but to interpret Fig. ^(a) prop- 
erly we must take into account the increase in both the 
dimension of configuration space and the average degree, 
(k), as the size of the network increases. 

For example, for a hypercubic lattice with a constant 
number of lattice points, L, along each edge and dimen- 
sion 3N, the number of lattice points, Matt, scales expo- 
nentially with N and the average number of steps be- 
tween any two lattice points is 3N{L + l)/3 = {L + 
1) log A'latt/S logL. By contrast, if the network behaves 
as a random graph, the average separation should scale as 

l0g^niin/l0g(fc) (X log Matt/ log(logiViatt) bcCauSC (fc) (X 

Ms/Mniin cx N. The sublogarithmic scaling suggested 
by Fig. ||(a) points to the latter scenario. This result 
is somewhat surprising. In Watts and Strogatz's small- 
world model the random-graph character results from the 
introduction of random links into the lattice, which can 




FIG. 3: The cumulative distribution for the number of nodes 
that have more than k connections. The curves correspond to 
clusters of different sizes, as labelled. An additional straight 
line with slope —(7 — 1), where 7=2.78, has been plotted to 
emphasise the power law tail. In the inset the cumulative 
probability distribution for the 12-, 13- and 14-atom clusters 
is plotted against k normalized by its average value, (k), to 
bring out the universal form of the distribution. 



potentially connect up distant nodes, but there is no ob- 
vious equivalent of the random links on the potential en- 
ergy landscapes. 

If we now examine the distributions for the numbers of 
connections for each node we find that as the size of the 
cluster increases a clear power-law tail develops, which 
has a universal form independent of the cluster size (Fig. 
^). The exponent of the power law, 7=2.78, is similar 
to other scale-free networks .Ej The cause of the random- 
graph like scaling of the average separation is thus the 
scale-free topology of the ISN. The network is extremely 
heterogeneous with a few hubs that have a very large 
number of connections, but with the majority of nodes 
only connected to a relatively small number of other min- 
ima. 

This is a particularly surprising result because all other 
scale-free networks are dynamic in origin. They grow and 
change over time, be it on an almost continuous basis 
as in the WWW or on evolutionary time scales in the 
case of biochemical networks. Even the recently intro- 
duced deterministic scale- fri ee . ri etworks are based on an 
iterative growth nrocedure.r^lr^l Furthermore, models of 
network growtljpp3 and studies on the time evolution of 
real networksBat2l strongly suggest that the heterogene- 
ity at the heart of the scale-free topology develops as a 
result of new nodes preferentially linking to those nodes 
which have many connections, be they much-cited papers 
or popular web-sites. However, the network associated 
with a potential energy landscape is static. It is simply 
determined by the form of the interatomic interactions 
and the number of atoms in the cluster. 

The source of the heterogeneity in the ISNs is apparent 
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potential energy of minimum 

FIG. 4: The dependence of the degree of a node on the po- 
tential energy of the corresponding minimum for the 14-atom 
cluster. The data points are for each individual minimum and 
the solid line is a binned average. 



from Fig. where we see that the degree of a node in- 
creases somewhat faster thaa exponentially as the energy 
of the minimum decreases-Ej The low-energy minima act 
as the hubs in the network. Thus, for the 14-atom clus- 
ter, 76% of the nodes in the network are connected to the 
global minimum. To measure the extent of the catchment 
basin around a minimum, we can calculate the distance 
in configuration space to all the transition states con- 
nected to a minimum. For LJ14 we find that this distance 
is 2.7 times larger for the global minimum than for the 
surrounding minima. When the multi-dimensionality of 
configuration space is taken into account this result sug- 
gests that the hyperarea of this catchment basin is many 
orders of magnitude larger than the average. Similarly, it 
has been previously found that on average the basin area 
falls off apppQximately exponentially with the energy of 
a minimum 

These results show that the global topology and to- 
pography of the potential energy landscape are inti- 
mately connected because the deep minima have very 
large catchment basins which are connected to lots of 
smaller basins around their edges. By contrast, if a po- 
tential energy landscape were flat and all basins of at- 
traction had the same area the scale-free topology would 
be lost. For example, an investigation of the network 
topology of thcj-configuration space of a non-interacting 
lattice polymeiEa found the connectivity distribution to 
be a Gaussian .Ej 

This contrasting example naturally leads one to ask 
how general is the topology that we have found for the 
Lennard- Jones clusters. Although such a question can 
only be definitively answered by similar analyses for a 
variety of systems, there is nothing "special" about the 
Lennard-Jones potential and so there is no reason why 
similar behaviour should not be seen for other materials. 



as long as there are no constraints present that would 
prevent the formation of the high degrees associated with 
the hubs. A polymer provides an example of the latter, 
because the connectivity of the chain limits the number 
of transition states that can surround a minimum. Fat 
example, a similar analysis for Lennard-Jones polymerstil 
did not find a power-law tail to the degree distribution 
(although it was still longer than exponential) O There 
is no equivalent of many of the transition states for the 
equivalent Lennard-Jones cluster because they involve 
the breaking of the polymer chain. 

The scale-free topology of the ISN is potentially good 
news for global optimization, the task of locating the 
global minimum. Even though the number of iniaima 
increases exponentially with the size of the system,E3 the 
average number of steps in the shortest path to the global 
minimum grows sublinearly with system size. Of course, 
finding this path is not necessarily easy. Our calculations 
of the shortest paths required information on the global 
structure of the potential energy landscape, whereas a 
global optimization algorithm usually takes a step based 
on only local information. 

Some path finding strategies to efficiently navigate 
scale-free networks have been suggested that make use 
of the fact that most of the-iiliprtest paths pass through 
the highly-connected hubsBS'LJ In our case the clear link 
between the topology and topography of the potential 
energy landscape provides an additional advantageous 
strategy. A downhill step to a lower-energy minimum is 
likely to take one to a minimum that is more connected 
and closer to the global minimum. How well-obeyed the 
latter correlation is, depends upon the global topography 
of the potential energy landscape and is a good indica- 
tor of the difficulty of global optimization. Thus, when 
the landscape is like a single funnel, global optimization 
algorithms can achieve near to the ideal scaling. For ex- 
ample, the basin-hopping algorithm can locate the global 
minimum of the 55-atom Lennard-Jones cluster after on 
average less than 150 minimizations when started from 
a random configuration, even though there are an esti- 
mated 10^^ minima.Eil The increasing number of links as 
the energy decreases evident in Fig. ^ further adds to 
the efficacy of a funnel in guiding the system towards its 
bottom, and provides clear evidence of the convergence 
of pathways into the hub at the funnel bottom that is of- 
ten postulated. By contrast, when an energy landscape 
has multiple funnels and there is a tendency to enter a 
funnel that leads the system to a low-energy minimum 
that is far from the global minimum, global optimization 
can be very difficult. 

The topology of the ISN will of course significantly af- 
fect the dynamics. This connection can be probed for 
very small systems where the network can be completely 
characterized and the inherent structure dynamics ob- 
tained by a master equation approach. However, this 
approach is not practical for the system sizes that are of 
most interest. Therefore, models of protein folding and 
the glass transition usually have to assume a simplified 
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topology for the interstate dynamics, or relate the 
dynamics to static quantities through phenomenological 
equations, such as the Adam-Gibbs equation which re- 
lates the relaxation tinae, in supercooled liquids to the 
configurational entropy.E2l To fully unlock the potential 
insights from the inherent structure view of the dynamics, 



a means of statistically modelling the network topology 
from a partial characterization of the potential energy 
landscape is therefore needed. Our results could signifi- 
cantly advance this goal. 
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